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Abstract 

A large calorimetric neutrino mass experiment using thermal detectors is ex- 
pected to play a crucial role in the challenge for directly assessing the neutrino 
mass. We discuss and compare here two approaches to the estimation of the ex- 
perimental sensitivity of such an experiment. The first method uses an analytic 
formulation and allows to readily obtain a sensible estimate over a wide range of 
experimental configurations. The second method is based on a frequentist Mon- 
tecarlo technique and is more precise and reliable. The Montecarlo approach is 
then exploited to study the main sources of systematic uncertainties peculiar to 
calorimetric experiments. Finally, the tools are applied to investigate the op- 
timal experimental configuration for a calorimetric experiment with Rhenium 
based thermal detectors. 
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1. Introduction 

Assessing the neutrino mass scale is one of the major challenges in today 
particle physics and astrophysics. This requires to measure the mass of one 
of the three neutrinos and the kinematical neutrino mass measurement is the 
only model independent method. In particular, the electron anti-neutrino mass 
can be measured by precisely analyzing the kinematics of electrons emitted in 
beta decays. In practice this means measuring the minimum energy carried 
away by the anti-neutrino, i.e. its rest mass, by observing the highest energy 
electrons emitted in the decay. To date, the study of the 3 H beta decay end- 
point by means of electrostatic spectrometers has proved to be the most sensitive 
approach, yielding an upper limit on the electron anti-neutrino mass of 2.2 eV 
Starting from 2012 the new experiment KATRIN will analyze the 3 H beta decay 
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end-point with a much more sensitive electrostatic spectrometer and with an 
expected statistical sensitivity of about 0.2 eV Q 

However, these spectrometric experiments suffer from many systematic un- 
certainties because the measured electron energy has to be corrected for the 
energy lost in exciting atomic and molecular states, in crossing the source, in 
scattering through the spectrometer, and more. To avoid these uncertainties it 
was proposed to embed the beta source in a detector and to perform a so called 
calorimetric measurement. Ideally, in such a configuration for each decay the 
detector measures all the energy released except for the energy carried away by 
the neutrino0 

A drawback of calorimetry is that a calorimeter is forced to detect all the 
beta decays while only the ones very close to the end-point E are useful for 
measuring the neutrino mass. The fraction of useful decays in a small interval 
AE below Eq is approximately given by (AE/Eq) 3 , therefore it pays out to 
select a beta decaying isotope with the lowest Eq value. In the past, calorimet- 
ric neutrino mass experiments have been performed implanting 3 H in Silicon 
diode detectors Q . The use of 187 Re as beta source seems more promising since 
it is the beta-active nuclide with the second lowest known transition energy 
(E Q ~ 2.5 keV). In the '80s S.Vitale proposed to realize calorimetric neutrino 
mass using 187 Re as beta source and exploiting the thermal detection tech- 
nique m. Since then two experiments of this kind has been carried out: the 
MANUyl 0] and MIBETAjl, @| experiments. MANU used one detector with 
a NTD thermistor glued to a 1.6 mg metallic rhenium single crystal, while MI- 
BETA used an array of ten silicon implanted thermistors with AgRe04 crystals 
for a total mass of about 2.2 mg. The two experiments collected statistics corre- 
sponding to 10 7 and 1.7x 10 7 decays respectively, yielding limits on m„ of about 
26 eV at 95% CL and 15 eV at 90% CL respectively. 

Recent developments in the thermal detection technique let think about a 
new very large calorimetric experiment aiming at a sub-eV sensitivity: this it 
what the MARE project is about [I, Ell- 
in this paper we present a comprehensive discussion of the potential sen- 
sitivity to the neutrino mass for a calorimetric experiment. First, through an 
analytical approach, we derive an algorithm to assess the statistical sensitivity 
for a given experimental configuration. Then a Montecarlo method is described 
which allows to get more precise statistical sensitivity estimates. The results of 
the analytic approach are then validated through the comparison with the Mon- 
tecarlo results over a wide range of experimental parameters. The second part 
of this paper focuses on the systematic uncertainties peculiar to the calorimetric 
technique by applying extensively the Montecarlo approach to their investiga- 
tion. We conclude with a discussion of the possible experimental configurations 
for future large scale calorimetric experiments. 



1 In practice particles emitted at the detector surface or with enough energy may also escape 
detection. In most cases, however, only a small fraction of the decays are affected by such an 
effect. 
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2. Statistical sensitivity 



2.1. The analytic approach 

In the following we derive an approximate analytic expression for the statis- 
tical sensitivity of a calorimetric neutrino mass experiment. The primary effect 
of a finite mass m v on the beta spectrum is to cause the spectrum to turn more 
sharply down to zero a distance m v below the end-point Eq (lower panel of Fig- 
ured]). To rule out a finite mass, we must be sensitive to the number of counts 
expected in this interval. The fraction of the total spectrum within an interval 
AE below the end-point Eq is given by 



E -AE 



N[)(E,m u )dE 



(1) 



where Np(E,m v ) is the beta energy spectrum for a neutrino mass m„ and 
normalized to unity. The signal to detect in counts is therefore 



signal = ApN det \F^{m v ) - Fae(0)^m 



(2) 



where Ap is the single detector source activity, A^et is the number of identical 
detectors and tu is the measuring time (see lower panel of Figured]). 
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Figure 1: Higher panel: beta spectrum as in JS) compared with pile-up spectrum (|3j|. Lower 
panel: zoom around the end point, with a comparison between and finite-neutrino-mass 
beta spectra 
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The detection of this signal is impaired by the noise caused by the statistical 
fluctuations of the total measured spectrum in the interval AE. For a calorime- 
ter the total measured spectrum is obtained summing up - from all detectors 
- the beta decay events, the counts due to unresolved pile-up of two or more 
decays, and any additional background counts. 

As a first approximation we can neglect the pile-up of more than two de- 
cays. We can then crudely approximate the pile-up spectrum by assuming a 
constant pulse-pair resolving time, tr, such that events with greater separation 
are always detected as being doubles, while those at smaller separations are 
always interpreted as singles with an apparent energy equal to the sum of the 
two events. In fact, the resolving time will depend on the amplitude of both 
events, and the sum amplitude will depend on the separation time and the filter 
used, so a proper calculation would have to be done as a Monte Carlo with 
the actual filters and pulse-pair detection algorithm being used. However, this 
approximation is good enough to get the correct scaling and an approximate 
answer. 

The parameter tr is related to the detector signal bandwith and high fre- 
quency signal-to-noise ratio: in practice tr is of the order of the detector rise 
time. 

The two event pile spectrum is given by 

Npp(E) = (1 - e- A ? TR )N p {E, 0) <g> Np(E, 0) (3) 

and the fraction of this spectrum within the interval AE below the end-point 
Eq is obtained by 

pEq pEq 

f ae = / N pp {E)dE a T R A fj / Np(E, 0) ® N P (E, 0)dE (4) 
JE -AE Je -ae 

where, at first order, trA/3 is the probability for the two event pile-up to occur, 
i.e. the fraction of unresolved pile-up events, f pp . From ([T|) and (01 one can 
write the noise in counts as 



noise 



yj ApNteiFusiO) + F™ E )t M + N de tbAEtM (5) 

where b is the average background counting rate for unit energy and for a single 
detector. We can then write the signal to noise ratio in a region within AE of 
the end-point Eo as 

noise V p y/F^)+F^ E + bAE/A p V ' 

It is now useful to introduce the exposure T = N,iettM & n d total number of 

events or total statistics of the experiment N ev = ApNdettM- 

The value of m„ which makes this ratio equal to 1.7 is the sensitivity at 90% 

confidence level, £90 (m„). Therefore one has to solve for m u the following 

equation 

rzr- \F±E{m v ) - f*e(q)\ =i7 



^FAE^ + F^ + bAE/Ap 
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To evaluate ([7]), we can consider approximate expressions for F^{m v ) and 
F^ E . In particular if we restrict ourselves to 187 Re, which has a first forbidden 
unique beta transition, we can make use of the following empirical spectrum 



3 / TTh^" 

ME,m v )^^(E -Er ] Jl- JE -^ w (8) 
which is an extremely good approximation of the expected theoretical shape 



11 1 as well as a perfect description of the experimental observations [1, 0] • 
For a null mass, from (151) we can derive 



AE N 3 



Fae(0) = (^) (9) 

while, for a small but finite mass m Ul using a second order expansion in m v / AE 
we have approximately 

^^^(i-l^ + lgr) do) 

For the pile-up spectrum ([3]) using ([8j we can calculate (between and E$) 

N„ (E) = (1 - e—)i- ( 9 | - 18 1 + 12 J - 3 J + A J) ,11) 

Between Eq and 2Eq the expression for N pp (E) is more complicated and it is of 
no use in this context. Substituting (fTT|) in ([4} and carrying out the integration, 
we obtain 



1 / AE AE 2 AE 3 AE 4 AE 5 AE 6 \ 
,— 6 h 15 — 5- + 20 — ^ - 15 — -r- - 6 — r «- 

20 V #0 El El El E% El J 



where we have used the approximation (1 — e A TR ) ss A^tr — f pp because in 
all interesting experimental configurations f pp <C 0.01. 

Substituting ©, dTU]) and (HU) in (0), keeping only the terms up to (AE/E ) 3 
and considering that (1 + f pp ) ~ 1, we obtain 



m2/3 Ar , 3ml\ r— / AE 3 , / 3A.E 3A£ 2 \ 

(13) 

which can be solved for m ll to give the sensitivity at 90% confidence level, 
£go(m„). Considering only the leading terms in (jTUJ) and (Tj"2"]) then the solution 
is 



S9o(m„) = 1.13 



N„ 



/\E Eq / 3 Eq 



E Q ' AE V10 + 



(14) 



In order to make meaningful use of (|14p one has to interpret correctly the energy 
interval AE. The same applies to the solution of (Q]: 
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The two terms in the square bracket in (|14[) represent the contributions 
to the noise from the statistical fluctuations of the beta and pile-up spectra 
respectively (here we neglect the background term for sake of clarity). When 
the pile- up term is negligible (because of the low rate Ap or of the short resolving 
time tr) the left term dominates and it pays out to keep AE as small as possible: 
the limit is of course the detector energy resolution. On the other extreme, when 
the end-point of the beta spectrum is buried in the pile-up spectrum, the noise 
is dominated by the right term. In this case the signal-to-noise ratio improves 
by enlarging the energy interval AE. 

It is then clear that there is no defined value of AE to plug in (fT?)) : the 
solution we have found is to choose AE as the value that minimizes T,go(m u ) 
for a given set of experimental parameters, with the boundary condition that 
AE cannot be smaller than the detector energy resolution A^fwhm- 

In particular, for the simpler case of (fl4|) . by searching the positive zero of 
the derivative with respect to AE, we can obtain the following 



This approach for defining AE can be applied to the solution of (fT3|) as well. 
In this case AE can be evaluated numerically following the above prescriptions. 
In particular to obtain the results presented in this paper, after making the 
substitution AE — > |A23opt| + AE'fwhm, we have found numerically the AE opt 
which minimizes the solution of (|13p 



It is worth noting that, in this analysis, equations (fT5|) and (p~6|) are the only 
places where the detector energy resolution AUfwhm shows up. 

2.2. Montecarlo approach 

In this section we describe a frequentist Montecarlo code developed to esti- 
mate the statistical sensitivity of a neutrino mass experiment performed with 
thermal calorimeters. The approach is to simulate the beta spectra that would 
be measured by a large number of experiments carried out in a given configura- 
tion: the spectra are then fit as the real ones @ and the statistical sensitivity 
is deduced from the distribution of the obtained parameters. 

The Montecarlo parameters describing the experimental configuration are 
the total statistics N ev , the FWHM of the Gaussian energy resolution ABfwhm, 
the fraction of unresolved pile-up events f pp and the background B(E). These 
input parameters can be derived from the ones actually characterizing a real 
experiment: N ev = N det Apt M and, by recalling (gj), f pp = F^ e=Eq « A t r , 
where again N^et is the number of detectors, Ap is the beta decay activity of a 
single detector, tu is the measuring time and tr is the pile-up resolving time. 

The procedure to estimate the statistical sensitivity goes through the follow- 
ing steps: 




(15) 



£go(»v) = / (AE opt , AEfwhm,tr, A/3, Ndet,tM,b) 



(16) 



G 



The theoretical spectrum S(E) which is expected to be measured by the 
virtual experiments is evaluated: 

S(E) = [N ev {N p {E,m v ) + f pp N (E,O) ® N (E,O)) + B{E)]® R(E) 

(17) 

where Np(E,m v ) is again the 187 Re beta spectrum normalized to unity 
(J5J), B(E) the background energy spectrum and R{E) is the detector en- 
ergy response function. The B(E) function is usually taken as a constant 
B(E) = bT. The response function R(E) is assumed to be a symmetric 
Gaussian 

G(E) = —^=e~^ (18) 



ctV2tt 

with standard deviation a = AEfwhm/2.35. 

The virtual outcome of a large number (between 100 and 1000) of exper- 
iments is numerically generated by letting the spectrum S(E) fluctuate 
according to a Poisson statistics. The simulated experimental spectra are 
generated on an energy interval which is smaller than the full - 2Eq 
interval. 

Each simulated spectrum is fitted using ()17|) and leaving m 2 , Eq, N ev , 
f pp and b as free paramaters. The fit is restricted to an energy interval 
smaller than the one used for the simulated spectrum generation. 

The 90% C.L. rriy statistical sensitivity £90(777^) of the simulated experi- 
mental configuration is given by £90 (m„) = ^/1.7<7 m 2 , where a m 2 is the 
standard deviation of the distribution of the m 2 found by fitting the spec- 
tra. 



l — J2« - ml) 2 = j^jW - ml*) (19) 



where N is the number of generated spectra and m 2 . are the values found 
in each fit for m 2 fit parameter. 

The statistical error on the 90% C.L. m u statistical sensitivity is estimated 
as follows. By defining yi = {m 2 . — m 2 ) 2 , we have y ~ and we can 
write 



2 #r5 2 N 

The error on y/lj = a m 2 is given by 

1 



9 \4 4 

mi) - al 



y Na m 2 

and therefore the error on £90(777,,,) is obtained 



(20) 



(21) 
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Figure 2: Comparison between the statistical sensitivity as estimated by a Montecarlo ap- 
proach (symbols) and by the analytic formulation (lines). The continuous and dashed lines 
are obtained using H14I I with H15I I and II161 respectively. The statistical sensitivity is evaluated 
for an exposure T = 10000 detector X year. The rightmost point corresponds to a total statis- 
tics N ev = 3 X 10 12 . The upper panel shows how AE op t gets larger as the pile-up spectrum 
increase its weight. 



Using equation one finds that the statistical error on the Montecarlo 
results is around 3% and 1% for about 100 and 1000 simulated experiments 
respectively. 

2.3. Analytic vs. Montecarlo 

We have compared the predictions of the two approaches described in the 
previous sections for a wide range of experimental configurations suitable for 
obtaining a sub-eV neutrino mass sensitivity (FigureEl [3j [4] and [5|) . In all plots 
the continuous lines are obtained through equations (fT4|) and (fT"5|) . while the 
dashed lines are obtained from (|T5)) as described at the end of the section on 
the statistical sensitivity. The symbols are the results of the Montecarlos, which 
have negligible errors on these scales (see equation (|22p). 

Figured shows how the sensitivity improves for increasing single detectors 
activity Ap, with the other experimental parameters fixed and for an exposure 
T of 10000 detx year. From this plot one can deduce that it pays out to increase 
the single detector activity at the expense of an increased fraction of pile-up 
events, f pp . The upper panel of Figure[2] shows the value of the optimal energy 
interval AE opt defined as discussed above. 
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Figure 3: Comparison between the statistical sensitivity as estimated by a Montecarlo ap- 
proach (symbols) and by the analytic formulation (lines) for a total statistics N ev of 10 14 
events. Continuous and dashed lines are as in Figure[2] 




Figure 4: Comparison between the statistical sensitivity as estimated by a Montecarlo ap- 
proach (symbols) and by the analytic formulation (lines). Continuous and dashed lines are as 
in Figure|2] The fine dashed line shows a functional dependence for sake of comparison. 
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integrated signal-to-background ratio 



Figure 5: Comparison between the statistical sensitivity as estimated by a Montecarlo ap- 
proach (symbols) and by the analytic formulation (lines). The integrated signal-to- background 
ratio is given by N ev /N),kg where Nhkg = bE$T (see text) and N ev = 10 14 . An integrated 
ratio of about 3 X 10 4 corresponds to the background level measured in the Milano experiment. 



Figure[3] demonstrates how the impact of the detector FWHM energy res- 
olution AE on the sensitivity is reduced by a high pile- up fraction f pp . This 
Figure shows also one of the limits of the analytic approach described above, 
i.e. the poor consideration of the detector energy resolution which translates in 
a too weak dependence of the sensitivity on this parameter. 

From Figure|4] it is possible to appreciate the importance of the total statis- 
tics N ev to reach a neutrino mass sensitivity of the order of 0.1 eV. In particular 
energy resolutions AE of about 1 eV and pile-up fractions of the order of 10 -6 
are required. The fine dashed line on the plot corresponds to a N e ^^ func- 
tional dependence of the sensitivity (recall equation (|14p ): this dependence may 
be exploited to scale the Montecarlo results. 

In Figurc[5] the impact of the continuous background below the beta spec- 
trum is shown. The abscissa is the ratio between N ev and the total number 
of background counts Nbk g between and E , i.e. Nbkg — bE T. Clearly the 
impact is lower for higher pile-up fractions f pp . 

For more details on each plots the reader can refer to their captions. In the 
final section of this paper these results will be used to assess the potential of a 
calorimetric neutrino mass experiment using 187 Re. 

Although the agreement is only partial, the comparison confirms that the 
analytic formulation goes in the right direction to be used to make useful predic- 
tions. Nevertheless we believe that the most accurate estimate of the sensitivity 
is the one obtained through the Montecarlo frequentist approach. From the 
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Figures it is apparent that the formula tends to overestimate the sensitivity for 
increasing AE opt , i.e. when the sensitivity is limited by the pile-up or the con- 
tinuous background. To improve the predictive power of the analytical approach 
it is possible to introduce free parameters in the formulas and to adjust them 
to better reproduce the Montecarlo predictions, but this is out of the scope of 
this paper. 

As a general check, both approaches can be applied taking as input the exper- 
imental parameters of the Milano experiment with an array of AgRe04 crystals 
whose results are presented in [f|. For N ev = 1.7 x 10 7 , A£?fwhm = 28.5eV, 
f pp = 2.3 x 10 -4 and N ev /Nf,k g — 3.28 x 10 4 one obtains a sensitivity at 90% C.L. 
of about 15 and 17 eV, using equations (fbl)) and the Montecarlo respectively, 
while the limit on the neutrino mass reported in [|| is about 16 eV at 90% C.L. 

3. Systematic uncertainties 

Although, as mentioned in the introduction, a calorimetric neutrino mass 
experiment is considered free from systematics related to the external source 
effects, still it may be affected by other uncertainties. 

As it will be discussed in the following, the origin of some of these effects 
(electron escape, beta decay spectral shape and beta environmental fine struc- 
ture) is indeed related to the beta source and may be unavoidable in spite of 
the calorimetric configuration. In order to minimize the related uncertainties, 
this kind of effects must be precisely modelled with the help of theoretical in- 
vestigations, independent experiments and Montecarlo simulations. 

Other systematic uncertainties arise from instrumental effects and can be 
mitigated through an improved detector design and characterization as well as 
a careful off-line data analysis. 

The frequentist Montecarlo code described above can be readily adapted to 
estimate the many systematic effects which in general fall under two categories. 
To the first belong the uncertainties due to lack of accuracy with which the ex- 
perimental parameters are determined. To assess the corresponding systematic 
uncertainties in the generated spectra, the parameters are randomly fluctuated 
- according to the given accuracy - while they are kept fixed to their average 
value in the fitting function S(E) (see f| 1 T|) ^) . The second category consists in 
the effects caused by an incomplete or incorrect modelling of the data. In this 
case the adopted approach is to include the effects in the generated spectra, but 
not in the fitting function S(E). 

In general these procedures result in a shift of rr? v away from zero and, in 
some cases, in a sensible deterioration of the sensitivity as shown by the wider 
error bars in the plots. The purpose of this analysis is to identify the size of the 
inaccuracy or of the neglected effect for which the shift of to 2 , remains within a 
given limit. TablefT] reports the results in terms of the systematic uncertainty 
giving a shift of |to 2 | less than 0.01 eV 2 - value for which the systematic uncer- 
tainties remain smaller than the statistical error of an experiment aiming at a 
sub-eV sensitivity. 
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In the following we analyze in more details the effects that are more likely to 
introduce systematic uncertainties in calorimetric neutrino mass experiments. 
Except where differently stated, the plots in this section are obtained for N ev = 
10 14 , A£: FW HM = 1.5 eV and f pp = 1CT 6 . 

3.1. Source related uncertainties 

Excited final states. Beta decays to any kind of excited final state are not going 
to perturb the beta spectrum end-point as long as the state lifetimes are shorter 
than the detector integration time - which is always more than about 1 fis. 
In fact, while an excitation energy E exc lost in the beta decay shifts the beta 
spectrum end-point to E' = Eq — E exc < Eq, at the same time the coincident 
detection of the energy E exc released in the state de-excitation adds an energy 
offset to the beta spectrum. The final outcome is that to each excited state 
corresponds a beta spectrum starting at E exc and with end-point Eq = E' Q + 
E exc . Thanks to its simple functional energy dependency ([8]) the 187 Re beta 
decay spectral shape above E exc remains unperturbed!! 

Electron escape. A fraction of electrons emitted in the decays of 187 Re nuclei 
close to the detector surface will not be contained in the calorimeter. We have 
used the Geant4 toolkit [3] to estimate the type and magnitude of this effect 
on the measured spectrum. Figure[6] shows the results for a 1 mg cubic Rhenium 
detector - i.e. with linear dimensions of about 0.362 mm - in terms of relative 
deviation with respect to the spectrum given by ([5]). The simulation has been 
repeated for the two available Geant4 low energy extensions (the results in 
Figure[J5] are the ones obtained using the Penelope extension) and for different 
low energy cuts applied in the electron transport. These tests has confirmed the 
shape of the effect and its magnitude, while giving slightly different results. In 
conclusion the Montecarlo simulation cannot be considered reliable to precisely 
calculate the effect, also considering the uncertainties in the shape and size of 
the detector rhenium absorber. The effect on the measured spectrum can be 
parametrized as a multiplicative factor to include in ((8|) given by 

f esc {E) = l-a esc ^- (23) 
Eo 

where the dimensionless a esc parameter will have to be left free in the data 
analysis. The solid line in Figure[f5] corresponds to a esc — 1.9 x 10~ 5 . We 
estimated the systematic error arising when this effect is not included in the 
data analysis for various values of a esc . The results are plotted in Figure[7l 



2 This is not true for a more general beta spectrum shape as for example in the case of 
allowed transitions with N(E) oc pE(Eo — E) 2 F(Z, E), where p is the electron momentum and 
F(Z, E) is the Fermi factor. For calorimeters, it is more generally true that at the end-point 
the correction for the presence of excited final states vanishes approximately as 5^(1 + Vi/E), 
where Vj are the final state energies. 
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electron energy [eV] 



Figure 6: Spectral distortion caused by surface electron escape as calculated by a Montccarlo 
simulation for a 1 mg cubic Rhenium detector. The effect plotted in the graph is (N'(E) — 
N(E))/N{E), where N(E) is given by equation © and N'(E) = N(E)f eB c(E) 
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Figure 7: Systematic shift caused by ignoring the correction for the surface electron escape. 
Points obtained for N ev = 10 14 , A_E FWHM = 1.5eV, f pp = 1CT 6 and bT = Oc/eV. 



3.2. Beta spectrum uncertainties 

Spectral shape. Although the use of equation ([5]) for the 187 Re beta decay spec- 
trum is up to the purpose of the present work, future high statistic experiments 
will need a more precise description of the spectrum. In order to estimate the 
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sensitivity to deviations from the simple equation ([5]) we have considered the 
corrective factor 

fcorr(E) = \ + ai E + a 2 E 2 (24) 

which is actually an extension of the correction for the escape of beta elec- 
trons (f2"3"| . The effect of deviations of this kind has been investigated by 
generating the experimental spectra according to a modified beta spectrum 
N'(E) = N(E)f corr (E) and fitting them using the regular beta spectrum N(E) 
in S(E) (fT7)l. Figure[5J shows the effect on m 2 varying the a, coefficients (see 
the caption for more details) . Figure[5] displays the results for positive values of 
the <Zj coefficients. 
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Figure 8: Systematic shift due to a deviation of the beta spectrum shape from the simple 
quadratic form in ([Sj - The four graphs represent the four quadrants of the (a\,a,2) plane, 
where ai and ai are the correction coefficients introduced in 11241 . Starting from upper right 
graph and going clockwise, they are the ai vs. 02, ai vs. — 02, — ai vs. — 02, and — a\ vs. 02 
plots. The colour coded z-axis is the value. The Montecarlo parameters are N ev = 10 14 , 



AEfwhm = 1.5 eV f p 



icr 



and bT = Oc/eV. 



For sake of comparison, one can consider that the deviation of the approx- 
imate beta spectrum shape described by equation {5J from the theoretical one 



given in llj can be parametrized as 



fcorr(E) ps 1.0 - 1.8 x 10- 5 E + 2.8 x l(T 10 i; 2 - 3.5 x 1CT 15 £; 3 + ... (25) 

Neglecting such a correction in the data analysis would systematically shift mf, 
by about -280 eV 2 . 

Beta Environmental Fine Structure. The Beta Environmental Fine Structure 
(BEFS) is a modulation of the beta emission probability due to the atomic and 
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Figure 9: Detail of the a\ > vs. 122 > quadrant for a wider parameter range and with 
logarithmic z-axis. 




electron energy [keV] 

Figure 10: Residuals from the fit of a Montecarlo generated beta spectrum with BEFS using a 
fit function without BEFS. The Montecarlo is for N ev = 10 10 , AB FWHM = 5eV, f pp = 10" 5 
and bT = c/eV. 



molecular surrounding of decaying nuclei: it is the analogous of the oscillation 
observed in the Extended X-ray Absorption Analysis (EXAFS) and it is ex- 
plained by the electron wave structure in terms of reflection and interference. 
Although the phenomenon is completely understood, its description is quite 
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complex and the parameters involved are not known a-priori. Because of its 
faintness, so far the BEFS has been observed in metallic Rhenium 12j and in 
AgRe04 [l3[ only below 1.5 keV where it is larger. It is clear that future neu- 
trino mass experiments will cope with the need of a very accurate description 
of the BEFS modulation up to the beta spectrum end-point. The parameters 
presently available are still affected by large statistical errors: for a safe extrap- 
olation up to the end-point the BEFS must be characterized using much higher 
statistics spectra. Meanwhile the Montecarlo approach can be used to show 



the shift on 



when data with BEFS included are fitted to a model without 



BEFS. For the BEFS function it is assumed that the one used to interpolate the 
data up to 1.5 keV can be used up to the end-point without modifications. A 
Montecarlo simulation of the Rhenium BEFS in a measurement with a statistics 
of about 10 10 events is shown in FigurefTUl in terms of residuals of the fit. Fig- 
ureQT] shows what happens to m} v when fitting spectra like the one in FigurefTOl 
with different left boundaries of the fitting energy interval. The effect worsens 
when the left boundary is moved to lower energies where the BEFS gets larger. 
The plot confirms that the inclusion of the BEFS in the end-point analysis is 
mandatory. 
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Figure 11: Systematic effect caused on by neglecting the BEFS when fitting the Montecarlo 
generated spectra with different left energy boundaries (see Figure JIO) . 



Pile-up spectrum. As discussed in § 12. 1[ formula §3§ holds only under the as- 
sumption of a constant resolving time tr. For real detectors tr depends on the 
pulse shape and on the noise level: in practice tr tends to increase for smaller 
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pulses. The detailed behavior of the resolving time is difficult to predict and 
must be modelled by Montecarlo methods taking in account both the actual 
pulse and noise frequency spectra and the algorithm used to identify the double 
pulses. Examples of such analysis can be found in 15, [l6j]. In order to get a 
sensible idea of the systematics related to this effect, we have used the results 
presented in [15j. We have numerically evaluated the pile-up spectrum intro- 



ducing a variable pile-up rejection efficiency described by an effective resolving 



time t 



eff 
R 



r-eff 



T R [l + (n T 



IV 



(26) 



where < 1 is the ratio between the amplitudes of the two events to discrim- 
inate, n TR is 2 and A r is 4.0. The function (|2"6"|) roughly approximates the one 
described in [l5j], even though it neglects the time ordering of the two events. 
With this approximation we find about 30% more unresolved pile-up events and 
a good description of the resulting pile-up spectrum for A^tr <C 0.1 is given by 



N'JE) = (1 - e- A '^)N p (E, 0) ® Np(E, 0)11 + 



0.35 



3 (B-E )/(480.0eV) _|_ i 



(27) 



FigurefT2l shows that the systematic shift caused by neglecting this deviation 
increases with the pile- up probability trA^. In particular it is apparent as a 
proper modelling of the pile-up spectrum is crucial for a pile-up rate as low as 
10~ 6 . 
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Figure 12: Systematic mjj shift caused by ignoring the pile-up spectrum correction for an 
energy dependent pile- up rejection efficiency. Points obtained for N ev = 10 14 , ABp^HM = 
1.5 eV, and bT = Oc/eV. 
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Figure 13: Instrumental systematic uncertainties for N ev = 10 14 , AEpwHM = 1-5 eV, f pp = 
10" 6 and bT = Oc/eV: response function uncertainty (left), energy calibration errors in an 
array (center) and detector energy resolution spread in an array (right). 



3.3. Instrumental uncertainties 

Response function uncertainty. In (|17j) the simplest response function R(E) 
used to model the data is a Gaussian G(E) (fTS|) which is completely determined 
by its standard deviation a = AE/2.35. The detector FWHM energy resolution 
AE is usually determined by means of a calibration procedure using radioactive 
sources. The accuracy with which AE, and therefore R(E), is known is mainly 
limited by statistics. 

Assuming a purely Gaussian R(E), the systematics due to the finite accuracy 
with which the detector FWHM energy resolution AE is known have been 
evaluated by letting fluctuate the detector energy resolution in the simulated 
spectra around a central value AE which is the fixed resolution used to fit 
the spectra. The detector energy resolution fluctuates according to a Gaussian 
distribution centered in AE with standard deviation a err (AE). The resulting 
shift of rv? v is shown in the left panel of Figure[T51 

With the number of events N ev and the pile-up fraction f pp considered in 
the present analysis, a calibration peak at an energy just above the beta decay 
end-point Eq would have only the pile-up spectrum as background. In these 
conditions, a perfectly Gaussian peak with 10 4 counts would allow an estimation 
of the FWHM energy resolution AE with an accuracy of about 1%. 

Response function tails. The actual response function R(E) may be as simple 
as a Gaussian, though presenting additional extra features which are difficult to 
identify in the calibration peaks. One example are small tails on the left side of 
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the main Gaussian peak. For a Gaussian with variance a, the function 



T(E) = A tai i — exp 



1 — erf 



E-En 



crA' 



aV2 V2; 



(28) 



represents an exponential tail with area A ta u and decay constant A. The effect 
of such an exponential tail has been studied using a response function R(E) = 
G(E) + T(E) in S(E) when generating the experimental spectra, while keeping 
the standard Gaussian response function in the fit: FigurefT4l shows the effect 
for various values of A and A ta ii < 1- 

Of course, identifying an exponential tail with a relative area A ta u as small 
as 10~ 4 requires a main Gaussian peak with statistics much larger than 10 4 
counts. 



-0.05 



> 

^ -0.1 



-0.15 



-0.2 



-0.25 



10' 





0.1 eV 1 


• •■>,- 


0.2 eV"' 


A-A X = 


0.5 eV"' 



1 



10 



10"" 



10"' 



tail area A 



tail 



Figure 14: Systematic effect caused on by an undetected exponential tail in a Gaussian 
response function R(E). The Montecarlo parameters are N ev = 10 14 , AEpwHM = 1.5eV, 



tp 



10" 6 and bT = Oc/eV. 



Calibration error in array. Future calorimetric experiments will be carried out 
with large arrays of thermal detectors (> 10 4 channels). Since each individual 
channel has to be energy calibrated prior to be summed up, a systematic effect 
may arise because of the finite calibration accuracy. For this analysis each simu- 
lated experimental spectrum is the sum of 10000 slightly mis-calibrated spectra. 
The mis-calibration is simulated shifting the energy calibration according to a 
Gaussian distribution centered in K with standard deviation <j(K), where K is 
the correct calibration factor. The center panel of Figurc[T51 displays the effect 
on ra 2 v . 

In the simplest hypothesis of a linear energy calibration without pedestal, 
a couple of calibration peaks close to beta end-point with about 10 4 counts 
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each and with a background arising solely from the pile-up spectrum would in 
principle allow to determine the calibration factor with an accuracy better than 



Response function dispersion in array. A second effect that can be observed 
when summing up many channels is due to the spread in the Gaussian response 
functions of the single detectors: the response function of the sum spectrum 
will not be Gaussian. In this case the simulated spectrum is the sum of 10000 
ones whose Gaussian response functions have FWHMs varying according to a 
Gaussian distribution centered in AE with standard deviation a sprea d{AE). 
The sum spectrum is analyzed assuming a response function with FWHM equal 
to AE. The results is shown in the left panel of Figure fT^I 

The dispersion of the energy resolution in an array strictly depends on the 
detector technology and it is therefore difficult to predict. Nevertheless, based 
on the experience with running arrays of thermal detectors, a o~ sprea d better 
than 10% should be realistic. 

Hidden background. In calorimetric experiments, since the beta source cannot 
be switched off, the background in the energy range of the beta spectrum cannot 
directly assessed. Therefore a costant background is usually included in the fit 
model S(E) as the safest hypothesis. Nevertheless we have analyzed the effect 
of neglecting this term. Figure[T5l shows the effect for various levels of constant 
background bT and confirms the importance of including the background term 



If no specific measure will be taken to reduce the background, future ex- 
periments are expected to have a signal-to-noise ratio similar to the one of the 
Milano experiment, i.e. about 10 4 -10 5 . In fact increasing the size of the single 
detector is unlikely to improve much the ratio since the background is expected 
to scale approximately as the detector mass. 

Eventually the background hidden below the beta spectrum could be not flat. 
We have explored this more critical situation making the simple hypothesis of 
a linear deviation from flatness starting just on the left of the pile-up spectrum 
end-point at 2E , expressed as 



The experimental spectra generated with the above linear background were fit 
with only the constant term in S(E): the results are shown in FigurefTOl for 
various values of constant background bT and various values for the deviation 
from flatness b\ . 



3 A similar effect due the spread in the resolving time tr is negligible since it affects only 
the total pile-up spectrum normalization. 



0.1%. 



in S(E). 




(29) 
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Figure 15: Systematic shift caused on by an undetected constant background bT under 
the beta spectrum. The Montccarlo parameters are N ev = 10 14 , AEfwhm = 1.5 cV and 
f pp = 10~ 6 . The Milano experiment signal-to-background ratio of about 3 X 10 4 would 
translates in to a constant background bT of about 10 6 c/eV on this plot. 
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Figure 16: Systematic shift caused on rrtj by an undetected deviation from flatness of 
the background under the beta spectrum. The Montecarlo parameters are N ev = 10 14 , 
A-EpwHM = 1-5 eV and f pp = 10~ 6 . The Milano experiment signal-to-background ratio cor- 
responds to a constant background bT of about 10 6 c/eV, therefore one order of magnitude 
higher than the worst constant background considered in this analysis. 
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Table 1: Analysis of systematic uncertainties for an experiment with N ev = 10 14 events, 
AE = 1.5 eV and f pp = 10 — 6 : for each uncertainty, the upper value which keeps the systematic 
shift of m 2 smaller than 0.01 eV 2 is given. 



source of 


quant ity 


maximum effect 


the uncertainty 


describing 


for 




the effect 


Ami < 0.01 eV 2 


electron escape 


Qesc 


1 x 10~ 5 


correction to 


K| (a 2 = 0) 


« 10- 9 eV~ 1 


theoretical spectral shape 


|a 2 | (ai = 0) 


» 10" 12 eV^ 2 


error on energy resolution AE 


a err (AE)/AE 


0.02 


tail in response function 
(A = 0.2 eV" 1 ) 


Afail 


1 x 10~ 4 


error on single channel energy 
calibration K 


<j{K)/K 


4 x 10" 4 


spread on energy resolution 
AE in the array 


a spread (AE)/AE 


0.1 


hidden costant background 


Nev/N bkg 


1 x 10 8 


hidden background linear 
deviation (bT = 10 5 c/eV) 


h 


w 0.1 



4. Future calorimetric experiments 

Given that the single channel activity A/3 is limited by technical considera- 
tions concerning the performance of the thermal detector (heat capacity, quasi- 
particle diffusion length, ...) the question is whether it is desirable to keep the 
pile-up negligible or not. There is no unique answer, although increasing the 
pile-up by increasing A/3 allows to accumulate more quickly large statistics, and, 
when pile-up dominates, the dependence on the energy resolution - which tends 
to degrade when A/3 increases - is attenuated (see Figured] and [3]). On the other 
hand the background caused at the end-point by the pile-up, together with a 
degraded energy resolution, may impair the ability to recognize and understand 
systematic effects. As a conclusion the optimal design of a neutrino mass ex- 
periment depends on the detection technique and, in particular, it depends on 
the effect of large absorbers on the detector performance. Nevertheless it may 
pay out to increase the single channel activity as much as possible, therefore 
relaxing the need of a high energy resolution. 

As an example, Tabled and report the scaled Montecarlo results for a 
target neutrino mass sensitivity £go(m„) equal to 0.2 eV and 0.1 eV respectively. 
Results are obtained in absence of background. The first line is a sort of baseline 
experimental configuration characterized by very demanding energy and time 
resolution and by very limited pile-up fraction f pp obtained by keeping the 
single detector activity A/3 at 1 Hz: in this conditions the target sensitivity is 
achieved with a relatively low statistics N ev at the expenses of a large required 
exposure T. In the other lines of the tables a larger activity of 10 Hz is considered 
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together with a progressive degradation of energy and time resolution. While a 
larger activity implies a lower required exposure, the poorer performances are 
compensated by the need of a larger statistics. From the tables it is clear that 
one can find a compromise between performances and exposure which is more 
convenient than the baseline high performance experimental configuration. 

For example a target neutrino mass sensitivity of 0.1 eV could be expected 
running for 10 years 3 x 10 5 rhenium detectors, each with a mass of 10 mg - 
giving an activity of about 10 Hz - and with energy and time resolutions of 
about 1 eV and 1 (is respectively. The total required mass of rhenium is about 
3 kg. 



Table 2: Experimental exposure required for a 0.2 eV m„ statistical sensitivity. 





TR 


AE 


N ev 


exposure T 


[Hz] 


H 


[eV] 


[counts] 


[detector x year] 


1 


i 


1 


0.2 x 10 14 


7.6 x 10 5 


10 


i 


1 


0.7 x 10 14 


2.1 x 10 b 


10 


3 


3 


1.3 x 10 14 


4.1 x 10 5 


10 


■5 


■5 


1.9 x 10 14 


6.1 x 10 5 


10 


10 


10 


3.3 x 10 14 


10.5 x 10 5 



Table 3: Experimental exposure required for a 0.1 eV m„ statistical sensitivity. 





TR 


AE 


N ev 


exposure T 


[Hz] 




[eV] 


[counts] 


[detector x year] 


1 


0.1 


0.1 


1.7 x 10 14 


5.4 x 10 B 


10 


0.1 


0.1 


5.3 x 10 14 


1.7 x 10 B 


10 


1 


1 


10.3 x 10 14 


3.3 x 10 6 


10 


3 


3 


21.4 x 10 14 


6.8 x 10 6 


10 


5 


5 


43.6 x 10 14 


13.9 x 10 6 



5. Conclusions 

In this paper we have thoroughly discussed the statistical sensitivity of calori- 
metric Rhenium based neutrino mass experiments. 

To estimate the statistical sensitivity, two methods have been developed. 
They are based respectively on an analytic and a Montecarlo approach: the 
results presented and compared in § 12.31 show the prominent importance of the 
total statistics collected by such an experiment in order to reach a sub-eV sen- 
sitivity. 
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Extending the application of the Montecarlo approach, we have then ana- 
lyzed the expected sources of systematic uncertainties peculiar to this kind of 
experiments. In particular, in jj l3.ll we have shown how crucial is for future ex- 
periments the understanding of the theoretical 187 Re beta decay spectrum and 
of BEFS. On the other hand, in ^ 13.31 we have shown how instrumental system- 
atic uncertainties may be kept under control by a proper characterization of the 
response function R{E) and by an accurate detector calibration: tasks that may 
be accomplished by controlling the calibration peak statistic. 

Finally, in §[U we have exploited the statistical analysis to devise a plausible 
experimental configuration capable to achieve a sensitivity of about 0.1 eV on 
the neutrino mass. 

As a concluding remarks, we believe we have demonstrated that calorimetric 
neutrino mass experiments with Rhenium based detectors offer a realistic chance 
to reach sensitivities comparable, or even beyond, the KATRIN goal. Moreover 
we have shown that, although systematics related to 187 Re beta decay theory 
and to BEFS require further investigations, these experiments should not be 
plagued by large systematic uncertainties. 

The authors wish to thank Prof. Dan McCammon for the many stimulating 
discussions on the topic of this paper. 

References 

[1] Ch. Kraus et al, Eur. Phys. J. C 40, 447 (2005). 

[2] KATRIN Desi gn Report (2004), FZKA7090; KATRIN Lol (2001), 
|hep-ex/0109033| 

[3] J.J.Simpson, Phys. Rev. D 23, 649 (1981). 

[4] S.Vitale et al., INFN Internal Report, INFN/BE-85/2, July 10, 1985. 
[5] C. Arnaboldi et al., Phys. Rev. Lett. 91, 161802 (2003). 
[6] M. Sisti et al., Nucl. Instrum. Methods A 520, 125 (2004). 
[7] F. Gatti, Nucl. Phys. B 91, 293 (2001). 
[8] M. Galeazzi et al., Phys. Rev. C 63, 014302 (2001). 
[9] the MARE proposal, |http:/ /mare. dfm.uninsubria.it 
[10] A. Nucciotti, J. Low Temp. Phys. 151, 597 (2008). 

[11] R. Dvornicky and F. Simkovic, "Measuring mass of neutrinos with beta 
decays of tritium and rhenium", talk presented at MEDEX'09, Prague, 
June 15-19, 2009. To be published in the proceedings. 

[12] F. Gatti et al., Nature 397, 137 (1999). 

[13] C. Arnaboldi et al., Phys. Rev. Lett. 96, 042503 (2006). 



24 



[14] https://geant4.web.cern.ch/geant4/ and the "Physics Reference Manual" 
therein. 

[15] F.Fontanelli et al., Nucl. Instrum. Methods A 421, 464 (1999). 
[16] L.Foggetta et al., J. Low Temp. Phys. 151, 613 (2008). 



25 



